The Effects of the Harmful Algal Bloom Species Karenia brevis on Survival of Red Porgy (Pagrus pagrus) Larvae

The harmful algal bloom species, Karenia brevis, forms annual, often intense blooms in the Gulf of Mexico, particularly along the west Florida shelf. Though the ability of K. brevis blooms to cause mass mortalities in juvenile fish are well documented, the direct effect of bloom concentrations on larval fish has not been studied extensively. To better understand the potential effect of K. brevis on larval fish survival, laboratory spawned red porgy (Pagrus pagrus) larvae from 4–26 days post-hatch were exposed to concentrations of K. brevis observed in the field for either 24 or 48 h. This species is representative of fish which spawn in regions of the Gulf of Mexico and whose larvae are epipelagic and may encounter K. brevis blooms. In this study, three different K. brevis strains varying in the amount of brevetoxin produced were tested. Larval survivorship was found to be inversely proportional to the amount of brevetoxin produced by each strain. The EC50 value from the combined 24 h experiments was ~163,000 K. brevis cells L−1, which corresponds to cell concentrations found in moderately dense blooms. Larval mortality also increased substantially in the 48 h versus 24 h exposure treatments. These findings indicate K. brevis blooms have the potential to contribute to natural mortality of fish larvae and further reduce inter-annual recruitment of fishery species whose stocks in the Gulf of Mexico may already be depleted.


Introduction
The harmful algal bloom species, Karenia brevis, forms annual and often intense blooms throughout the Gulf of Mexico [1]. This toxic dinoflagellate species produces brevetoxins (BTx) which are potent lipophilic polyether neurotoxins that disrupt nerve function in vertebrates by activating voltage-gated sodium channels [2,3]. Though the majority (~90%) of the BTxs are retained within the cells during a bloom [4], they can nonetheless rapidly permeate the food chain. This occurs primarily when filter-feeding bivalves,

Effects of K. brevis on Larval Fish Survival
The data from individual triplicate exposure experiments of 6 to 21 dph larvae exposed to equivalent cell concentrations of either strain NCMA2228 or SP3 were also compared (Table S1). There was no evidence for any difference in the susceptibility of earlier versus later dph (non-yolk sac) larvae to the toxic effects of K. brevis. Consequently, the data for the different dph larvae exposed to a specific K. brevis strain in a given year were combined. The mean survival of red porgy larvae was found to be inversely proportional to BTx concentrations per cell in the three different K. brevis isolates tested (Figure 1). Exposure to SP1, which contains only trace amounts of BTx, caused an initial decline in survival that was the same regardless of the number of K. brevis. After the initial uniform drop, survivorship remained stable and did not decline as K. brevis cell concentrations increased ( Figure 1A). NCMA2228 had an intermediate toxicity averaging 9.9 ± 2.2 pg BTx-1, -2, -3 cell −1 (n = 6) in the 2011 experiments and 10.9 ± 1.5 pg BTx-1, -2, -3 cell −1 (n = 3) in 2012 (Table S1). In both years survivorship declined in proportion to the number of cells present. The decline in survivorship versus cell concentration in 2012 was more pronounced when the BTx concentrations per cell were slightly higher ( Figure 1B). The highest toxicity strain SP3, averaged 14.5 ± 1.3 pg BTx-1, -2, -3 cell −1 (n = 2) and caused an even more rapid decline in survivorship than exposure to either SP1 or NCMA2228 for an equivalent number of cells ( Figure 1C). The only inter-year difference observed was the greater mortality of larvae in 2012 exposed to NCMA2228 compared to those exposed in 2011. A two-tailed t-test showed the difference between NCMA228 and SP3 at 100,000 cells for all replicates was not significant in 2012 (t (16) = 0.18, ns), but highly significant in the 2011 experiments (t (10) = 8.88, p < 0.001).

EC50 of Karenia brevis Versus Survivorship
The estimated EC50 value was 163,000 K. brevis cells L −1 (Figure 2). The average toxicity for the K. brevis isolates used to calculate this EC50 value was 9.4 ± 4.8 pg BTx cell −1 .

EC 50 of Karenia brevis Versus Survivorship
The estimated EC 50 value was 163,000 K. brevis cells L −1 (Figure 2). The average toxicity for the K. brevis isolates used to calculate this EC 50 value was 9.4 ± 4.8 pg BTx cell −1 .
Toxins 2022, 14, x FOR PEER REVIEW 5 of 15 Figure 2. Dose-response curve created using combined data for all the 24 h experiments where red porgy larvae were exposed to varying cell concentrations of the three K. brevis strains tested. The black dots represent the mean survivorship for the different days post hatch larvae exposed to a given concentration of Karenia brevis ± 1 SD. The resulting estimated EC50 for the red porgy larvae was ~163,000 K. brevis cells L −1 .
The corresponding estimated 24 h EC50 value based on the average BTx concentrations present in the different concentrations of whole cells tested was 17.8 pg L −1 ( Figure  S1). The data from the field indicate about 10% of the internal toxin pool is actually released into the water surrounding the cells [4].

Effect of Increased Exposure Time on Survivorship
Three trials were conducted where survivorship of the larvae exposed to NCMA2228 K. brevis were measured again at 48 h ( Figure 3) to determine the extent to which increased exposure time caused elevated mortalities. The dph selected are representative of range in larval ages tested and include from 4 dph larvae which still retained their yolk sac and two later post yolk sac ages. The results consistently showed a significant decline between 24 and 48 h exposure times in the 4 dph (97% to 76%: t(22) = 2.94, p = 0.008), 6 dph (85% to 47%: t(22) = 3.72, p = 0.001), and 16 dph (e.g., 73% to 40%: t(22) = 2.75, p = 0.011) larvae. Dose-response curve created using combined data for all the 24 h experiments where red porgy larvae were exposed to varying cell concentrations of the three K. brevis strains tested. The black dots represent the mean survivorship for the different days post hatch larvae exposed to a given concentration of Karenia brevis ± 1 SD. The resulting estimated EC 50 for the red porgy larvae was~163,000 K. brevis cells L −1 .
The corresponding estimated 24 h EC 50 value based on the average BTx concentrations present in the different concentrations of whole cells tested was 17.8 pg L −1 ( Figure S1). The data from the field indicate about 10% of the internal toxin pool is actually released into the water surrounding the cells [4].

Effect of Increased Exposure Time on Survivorship
Three trials were conducted where survivorship of the larvae exposed to NCMA2228 K. brevis were measured again at 48 h ( Figure 3) to determine the extent to which increased exposure time caused elevated mortalities. The dph selected are representative of range in larval ages tested and include from 4 dph larvae which still retained their yolk sac and two later post yolk sac ages. The results consistently showed a significant decline between 24 and 48 h exposure times in the 4 dph (97% to 76%: t (22) = 2.94, p = 0.008), 6 dph (85% to 47%: t (22) = 3.72, p = 0.001), and 16 dph (e.g., 73% to 40%: t (22) = 2.75, p = 0.011) larvae.

Discussion
Karenia brevis cell densities equivalent to those found in moderate to high density blooms were found to be highly toxic to red porgy larvae (Pagrus pagrus) (Figures 1 and  2). Red porgy is a broadcast spawner which is widely distributed over the west Florida

Discussion
Karenia brevis cell densities equivalent to those found in moderate to high density blooms were found to be highly toxic to red porgy larvae (Pagrus pagrus) (Figures 1 and 2). Red porgy is a broadcast spawner which is widely distributed over the west Florida shelf, most commonly at depths between 10 to 80 m [33,34]. Larvae tend to migrate into the top 5 m of the water column where they can readily come into contact with K. brevis blooms. This species also exhibits a high degree of site fidelity and spawns from December through March. [35]. K. brevis blooms usually initiate in the late summer or early fall and typically persist over significant portions of the mid to inner west Florida shelf until mid to late-winter (January-February) [1,36]. Though less common, intense K. brevis blooms sometimes fail to decline in late winter and instead continue unabated through the spring and following summer [1]. K. brevis cells are highly mobile, often migrating up in the water column during the day and down at night. When upwelling conditions are favorable, K. brevis cells can also rapidly accumulate hydrodynamically over a few days to a week reaching concentrations far exceeding the EC 50 values determined in this study [36,37]. Given their overleaping temporal and spatial distribution patterns, red porgy represent a good model species for evaluating the potential adverse impacts of K. brevis blooms on larval fish in the Gulf of Mexico.
The EC 50 value of 163,000 K. brevis cells L −1 obtained in this study indicates the greatest red porgy larvae mortalities will occur at moderate to high bloom densities (100,000-10,000,000 cells L −1 ; Figure 2). This estimate is consistent with K. brevis cell densities (≥100,000 cells L −1 ) known to predictably kill diverse juvenile and adult fish species [1,30,38,39]. The corresponding EC 50 value based on average BTx concentration in the cells was 17.8 pg L −1 . Interestingly, shellfish larvae appear to require higher cellular concentrations of Karenia, in the range of 1,000,000 to 3,000,000 cells L −1 , to induce significant mortality [40,41]. Though the EC 50 value obtained in this study was markedly low compared to those obtained for shellfish larvae exposed to K. brevis, it may still represent an overestimate. This is because the cellular BTx levels for the three K. brevis strains used in this study averaged 9.4 ± 4.8 pg BTx cell −1 compared to the measured mean toxicity of K. brevis blooms in field of 18.2 ± 5.5 pg BTx cell −1 [4]. The cause of these toxicity difference is primarily due to our experiments being conducted using exponentially growing cells. These fast-growing, non-nutrient limited cells exhibit lower toxicities than the slower growing, nutrient limited cells often encountered in the field [4,42]. For example, NCMA2228 averages 9 pg BTx-1, -2, -3 cell −1 in exponential phase growth, 19.2 pg BTx-1, -2, -3 cell −1 under nitrate limitation and 42.7 pg BTx-1,2,3 cell −1 under phosphate limitation [43,44]. Similarly, SP3 averages 13 pg BTx-1, -2, -3 cell-1 when in exponential nutrient replete growth and 19.7 when nitrogen limited [43]. Consequently, had field K. brevis populations containing almost two-fold higher cellular BTx concentrations been used instead of the strains employed in this study, a lower EC 50 estimate would be predicted. If true, exposure to even relatively low-density K. brevis blooms (50,000 cells L −1 or less) may cause increased larval mortalities.
Results from this study also showed that as exposure time was increased from 24 to 48 h, red porgy percent mortalities increased substantially regardless of larval age (Figure 3). This indicates the longer the larvae are exposed to relatively low K. brevis cell concentrations in the field, the greater the resultant larval mortalities. Riley et al. [31] noted 21 dph red drum larvae needed longer exposure times to reach the same level of mortality compared to yolk-sac larvae. The data on K. brevis cell concentration versus mortality data for the 6-21 post yolk-sac larvae indicated each age was equally susceptible to the toxic effects of K. brevis (Table S1). The data for the 4 dph yolk-sac larvae, in contrast, were insufficient to determine if these larvae were more or less susceptible to K. brevis toxicity compared to the post yolk-sac larvae.
Given the high toxicity of K. brevis and its wide distribution along the west Florida shelf, it is likely this toxic dinoflagellate may periodically cause larval mortalities sufficient to reduce inter-annual recruitment and year-class strength of important fish stocks. Any adverse population and community level associated effects would be compounded by the concomitant increases in juvenile and adult fish mortalities associated with K. brevis blooms [30,[45][46][47]. Further support for the potential negative impacts on fish stocks comes from the reduction in bivalve larvae recruitment [39,[48][49][50] and mass mortality of benthic organisms including corals observed during bloom during K. brevis blooms [22,45,[51][52][53][54]. In the latter case, radical shifts in the benthic community structure persisting for over two years following a major bloom have been observed. Any significant reductions in year-class strength are of particular concern for already overfished species in the Gulf of Mexico. Cumulatively, these data suggest more detailed correlative studies examining decadal changes in the nekton community and year-class strength of fish species following moderate to severe K. brevis bloom years are needed to fully elucidate the effects of bloom on the resilience of fisheries stocks in this region [46,55].
Though the toxicity of the K. brevis isolates correlated with their ability to produce different levels of BTx, it is important to note no experiments were conducted using purified BTx to conclusively demonstrate these neurotoxins per se were definitively responsible for the observed mortalities. Early histopathological studies indicated many fish exposed to K. brevis died from chronic tissue damage rather than neuro-intoxication [38]. One of the most likely candidates for such damage would be the production of reactive oxygen species (ROS) superoxide and hydrogen peroxide, which is known to be the mechanism of toxicity for some harmful algal bloom species [56]. Karenia brevis cells are known to produce low levels of intracellular ROS which are reported to both increase in concentration with culture age and fluctuated in a phasing diel pattern [57]. Other work confirmed K. brevis does indeed produce intermediate levels of extracellular ROS, but in contrast showed the production varied inversely with cell density [58]. The later study also failed to detect active ROS in spent media of K. brevis cultures. Recent surveys, while documenting the importance of ROS toxicity in K. mikimotoi, did not identify ROS as a major mechanism by which K. brevis causes toxicity [56,59]. Cumulatively, these data indicate the ROS is unlikely to account for the elevated larval mortality observed with increasing K. brevis concentrations.
Instead, K. brevis toxicity is widely attributed to BTXs and their derivatives being neurotoxic, hemolytic, and potent inducers of oxidative stress [60][61][62][63][64][65]. Its specific mode of action in blocking voltage gated sodium channels of most organisms in particular [28,66], suggests BTxs as the most parsimonious explanation for the observed toxicity of K. brevis to larval red porgy [60]. Regardless of the specific toxins involved, the results of this study indicate moderate to intense blooms of K. brevis may cause larval fish mortalities sufficient to affect population dynamics of ecologically and commercially important fish stocks in the Gulf of Mexico.

Larval Fish Rearing
Experiments were conducted using red porgy (Pagrus pagrus) larvae reared from two volitional spawns of brood stock held at the NOAA Beaufort Laboratory, Beaufort, NC, USA; one in 2011 and one in 2012. Viable eggs were separated by placing them in a 1 L graduated cylinder and pouring off the viable eggs into a 100 L hatching cone kept at 19 • C by a 3 4 horsepower heater/chiller. The larvae hatched 48 h later and were concentrated into a smaller cone via a siphon. Larvae were placed in 100 L recirculating (400% daily exchange) tanks at 23 • C and salinity of 34 [67]. Photoperiod was maintained at 16:8 light:dark cycle (07:00-23:00) from a single 90-watt Sylvania daylight spectrum fluorescent bulb. Aeration was supplied via a single flexible air diffuser around the center of a 500 µm mesh screen overflow. Larvae were fed S-type rotifers (Brachionus rotundiformis) enriched with Rotigrow plus (Reed Mariculture; Campbell, CA USA) starting 2 days post hatch (dph) at a density of 3 rotifers mL −1 and increased daily until 12 rotifers mL −1 was reached. Feeding took place once per day (08:00) by hand and an additional three times per day (12:00, 16:00, and 20:00) with refrigerated (10 • C) feed pumped into the tanks via a peristaltic pump controlled by a timer. A non-viable, condensed algae paste (Nanno, Reed Mariculture; Campbell, CA USA) was added twice daily at a concentration of 750,000 cells mL −1 starting at 1 dph as background. Instar I stage brine shrimp (Artemia salina) were co-fed with rotifers starting 15 dph and until 20 dph when Instar II stage enriched Artemia (DC DHA SeECo, INVE, Salt Lake City, UT, USA) were fed exclusively. Artemia were enriched 24 h prior to feeding as prescribed by the manufacturer's protocols. Artemia (1-3 Artemia mL −1 ) were fed the same way as the rotifers. Otohime C1 Fish Diet (Reed Mariculture; Campbell, CA, USA) was fed starting at 15 dph in addition to live organisms.
Red porgy (Sparidae; Pagrus pagrus) was selected as a representative species that exhibits an epipelagic larval stage with high likelihood of natural exposure to BTx in coastal marine waters [35]. Red porgy exhibit a two-phase life cycle where more sedentary and demersal adults produce pelagic larvae that are widely dispersive by hydrodynamic forces such as wind, tide, and turbulence [67]. By comparing larvae hatched in 2011 versus 2012, it was possible to evaluate if inter-year differences in the reared larvae affected susceptibility to K. brevis toxicity. The range of K. brevis concentrations tested represent bloom density levels from background to moderate as defined by the Florida Fish and Wildlife Research Institute [1]. The larval ages tested coincided with key developmental stages including the development of mouthparts and gills.

Karenia brevis Cultures
Three K. brevis strains, NCMA2228, SP3, and SP1, were used in this study. NCMA2228 was obtained from the National Center for Marine Algae and Microbiota and was isolated from Mote Marine Laboratory New Pass Dock, Sarasota, FL, USA. Strains SP1 and SP3 were provided by Dr. Edward Buskey of the University of Texas Marine Science Institute, Port Aransas TX, USA. Both SP1 and SP3 were isolated during a bloom event off the coast of Brownsville, Texas in October 1999 by S. Pargee, University of Texas Marine Science Institute, Port Aransas. K. brevis isolates have long been known to vary in toxicity [68] and these strains were selected because they are known to contain either trace (SP1), intermediate (NCMA2228), or high (SP3) intercellular levels of BTx [43,68,69].
Karenia brevis cells were grown in 2.5 L semi-continuous batch cultures as described in Hardison et al. [41]. . Nutrient-replete culture media contained 64 µM NaNO 3 and 4 µM NaH 2 PO 4 (N:P = 16:1, the Redfield ratio). Media were sterilized by microwave treatment [41]. Culture pH was measured initially and throughout each experiment with a Thermo Scientific (Waltham, MA, USA) Orion 3 Star pH meter equipped with a Ross ultra-combination pH electrode to ensure no carbon dioxide limitation occurred. Culture pH ranged from 8.10 to 8.30.
Cell concentrations were measured with a Beckman Coulter Inc. (Brea, CA, USA) Multisizer 3 electronic particle counter (0.5-mL sample volume) equipped with a 100-µm aperture. Cell growth curves were constructed as semi-log plots of total biovolume (µLcells Lmedia −1 ) versus time in days. Total biovolume (µLcells Lmedia −1 ) was calculated by multiplying the cellular concentration (cells L −1 ) by the mean volume per cell (µm 3 cell −1 ). Specific growth rates were computed from linear regressions of the natural log of total cell biovolume versus time after correcting for culture dilution [70]. Only exponentially growing, nutrient sufficient K. brevis cells were harvested and used in the exposure studies.

Brevetoxin Extractions and Quantification
Each time cells were used in an exposure experiment, aliquots of cell cultures were mixed 1:1 by volume with ethyl acetate for brevetoxin (BTx) determination using the methods described in Hardison et al. [43]. The ethyl acetate mixtures were sonicated for three minutes with a Branson Sonifier 250 (Branson Ultrasonics Corporation, Danbury, CT, USA) and complete cell disruption was confirmed by microscopy. The aliquots were desalted with Milli-Q water (Millipore Sigma, Burlington, MA, USA) and concentrated with a rotovap. Extraction efficiency was determined in each aliquot by the addition of an internal standard BTx-42-acetate (kindly provided by Dr. Dan Baden, University of North Carolina at Wilmington, Wilmington, NC, USA) and typically ranged from 95-100% (Lekan and Tomas, 2010). Brevetoxin concentrations in each aliquot were measured using an Agilent 1100 LC (Agilent Technologies, Santa Clara, CA, USA) coupled to a Thermo Finnigan TSQ Quantum triple quadrupole mass spectrometer (Thermo Fisher Scientific, Waltham, MA, USA), with an electrospray ion source interface. The LC-MS-MS methodology used has been described in detail [71,72]. An external standard curve of three purified brevetoxins (BTx-1, -2, -3) (World Ocean Solutions, Wilmington, NC, USA) was used to quantify the amounts of BTx-1, -2, -3 in each of the aliquots for each trial. Brevetoxin levels in this study are the mean total brevetoxin levels (pg cell −1 ), which were determined by taking an average of the sum BTx-1, -2, -3 across six aliquots in each trial. The main reason for measuring toxicity each time cells were used is that toxicity can vary even when grown under constant conditions [64].

Exposure Experiments
Experimental systems consisted of 1 L Pyrex beakers placed on the laboratory bench at room temperature (21 • C), full strength filtered seawater (35), and 24 h overhead fluorescent lighting. Experimental controls contained filtered seawater, larvae, and age-appropriate food, while experimental treatments contained filtered seawater, larvae, age-appropriate food source, and K. brevis. Three beakers (triplicates) of seven to ten larvae were used in each control and treatment. In 2011, larvae were initially tested at 4 dph using K. brevis strain NCMA2228 at concentrations from 0-115,000 cells L −1 ; however, survival was high across treatments, so a higher range of cell concentrations was tested in the remaining trials. Subsequently, larvae were exposed to NCMA2228 at 6,8,14,16, and 21 dph at K. brevis cell concentrations from 0-500,000 or 0-1,000,000 cells L −1 depending on the toxicity of the K. brevis strain being employed. Similarly, 13 dph larvae were exposed to K. brevis strain SP3 at concentrations from 0-1,000,000 cells L −1 . In 2012, preliminary trials with 5 dph larvae indicated the larvae were more sensitive to the slightly more toxic K. brevis strains used, when compared to the 2011 experiments. In this case, a lower range of cell concentrations were tested. Specifically, 7, 11, and 26 dph were exposed to K. brevis strain NCMA2228, 11 dph larvae to strain SP3, and 18 and 21 dph to strain SP1 at concentrations ranging from 0-500,000 cells L −1 . Cell concentrations were maintained throughout the duration of each trial and care was taken to prevent K. brevis cell damage and inadvertent release of excess toxin, which was confirmed by microscopically examining the cells to ensure they were healthy and intact. The average cell concentration for each of the three K. brevis isolates were then plotted versus percent survival ± one standard deviation and the toxicity of the isolate for each experiment noted.

Estimating the EC 50 Value
The 24 h-EC 50 value for larval red porgy exposed to different concentrations of K. brevis was calculated using the Graphpad Prism ® (GraphPad Software, Inc., La Jolla, CA, USA) program. Input data were the K. brevis cell numbers and the corresponding average 24 h survivorship ± one standard deviation for all of the experiments. Support for combining the data from multiple isolates having different BTx levels to estimate an average EC 50 value comes from a field study by Tester et al. [4]. They found an average toxicity of 18.2 ± 5.5 pg BTx per K. brevis cell −1 during a bloom event off the west coast of Florida. The measured toxicity of the K. brevis isolates used for each experimental treatment in this study was averaged to determine how closely their mean toxicity matched BTx concentrations per cell measured in the field. It should be noted that in addition to simulating average field toxicity, combining larval survivorship data for strains producing different amounts of BTX also increased the observed variance in percent survivorship for each K. brevis cell concentration tested.
An 24 h-EC 50 value for larval red porgy exposure based on pg BTx L −1 was also calculated as above using the same percent survivorship data. The only difference was that mean toxin content for each of the cell concentrations variously tested in the 24-h studies (0, 10,000, 100,000, 250,000, 500,000, or 1,000,000 cells L −1 ) was used instead of cell concertation as input to the Graphpad Prism program. The BTX values used to calculate the mean toxin values for each concentration of K. brevis cells tested were obtained by multiplying the measured BTx concentration (pg cell −1 ) for the strain used in each individual experiment by the specific concentration of Karenia cells L −1 being tested. The BTx concentrations were then sorted into groups corresponding to the cell concentration tested and corresponding average BTx concentrations calculated.

Examination of How Prolonged Exposure to K. brevis Impacts Survivorship of Red Porgy Larvae
For a subset of the experiments using K. brevis strain NCMA2228, the incubations with 4, 6, and 16 dph larvae from 2011 were allowed to proceed for 48 h and the resulting mortalities recorded (Figure 1). The toxicity of the isolates employed in these experiments were similar ranging from 11 to 12.3 pg BTx cell −1 . This provided a consistency in K. brevis concentration that allowed an estimate of how survivorship might be affected by exposure to toxic K. brevis blooms over longer periods of time. Statistical differences in mean survivorship between the 24 and 48 h exposure treatments was determined using a t-test.

Supplementary Materials:
The following supporting information can be downloaded at: https:// www.mdpi.com/article/10.3390/toxins14070439/s1, Table S1, Summary data for each of the 24 h and 48 h experiments carried out in this study where red porgy larvae were variously exposed to three different Karenia brevis isolates. Figure S1: Dose-response curve created using combined data for all the 24 h experiments where red porgy larvae were exposed to concentrations of the three K. brevis strains which contained different internal BTx concentrations (Table S1). The average BTx concentration for each concentration of K. brevis cells testes were calculated and used in conjunction with their corresponding mean survival rates (±SD) to estimate a BTx-based EC50. The resulting EC50 value was~17.8 pg L -1 in whole K. brevis cells. Field studies indicated about 10% of the internal BTx concentrations are released into the water column [4].  Institutional Review Board Statement: The animal study protocol was approved by Centers for Coastal Ocean Science for studies involving animals which follows the Federal Ethics Committee on Non-Human Biotechnology. Those guidelines deem acceptable research on fish whose results may be used to either benefit the fish themselves or to gain greater profit from fish thanks to the knowledge gained.

Informed Consent Statement: Not applicable.
Data Availability Statement: Not applicable. All data are shown in the manuscript or in supplementary Table S1.

Conflicts of Interest:
The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript; or in the decision to publish the results.